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Abstract 

We study the qualitative properties of cosmological models in scalar- 
tensor theories of gravity by exploiting the formal equivalence of these 
theories with general relativity minimally coupled to a scalar field 
under a conformal transformation and field redefinition. In particu- 
lar, we investigate the asymptotic behaviour of spatially homogeneous 
cosmological models in a class of scalar-tensor theories which are con- 
formally equivalent to general relativistic Bianchi cosmologies with 
a scalar field and an exponential potential whose qualitative features 
have been studied previously. Particular attention is focussed on those 
scalar-tensor theory cosmological models, which are shown to be self- 
similar, that correspond to general relativistic models that play an 
important role in describing the asymptotic behaviour of more gen- 
eral models (e.g., those cosmological models that act as early-time and 
late-time attractors). 



1 Introduction 



Scalar-tensor theories of gravity are currently of great interest, partially 
due to the fact that such theories occur as the low-energy limit in superstring 
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theory (see ]T| and references therein). The first scalar-tensor theories to 
appear (with u = u ) were due to Jordan [|, [J, Fierz and Brans and Dicke 
H and the most general scalar-tensor theories were formulated by Bergmann 
0, Nordtvedt and Wagoner ||. The observational limits on scalar-tensor 
theories include solar system tests || 0, |12| and cosmological tests such 
as Big Bang nucleosynthesis constraints [I 



m. 



The possible isotropization of spatially homogeneous cosmological models 
in scalar-tensor theories has been studied previously. For example, Chauvet 



and Cervantes-Cota |L5| have studied the possible isotropization of Bianchi 
models of types I, V and IX within the context of Brans-Dicke theory with- 
out a scalar potential, but with barotropic matter, p = (7 — by studying 
exact solutions at late times. Mimoso and Wands \TE\ have studied Brans- 
Dicke theory with a variable to = u(<f>) in the presence of barotropic matter 
(but with no scalar field potential) and, in particular, gave forms for ui under 
which Bianchi / models isotropize. We note that there is a formal equiva- 
lence between such a theory (with 7 7^ 2) and a scalar-tensor theory with a 
potential but without matter, via the field redefinitions V = (2 — 7)/^ and 

In a recent paper (see also |L8| and |0|), cosmological models contain- 
ing a scalar field with an exponential potential were studied. In particular, 
the asymptotic properties of the spatially homogeneous Bianchi models, and 
especially their possible isotropization and inflation, were investigated. Part 
of the motivation for studying such models is that they can arise naturally 

for example, 



in alternative theories of gravity [20 



21] has shown 



Halliwell 

that the dimensional reduction of higher- dimensional cosmologies leads to an 
effective four-dimensional theory coupled to a scalar field with an exponential 
self-interacting potential. 

The action for the general class of scalar-tensor theories (in the so-called 
Jordan frame) is given by || |8|, 



S 



-9 



'^R-^9 a \ 



2V{<f>) + 2C r 



d x. 



However, under the conformal transformation and field redefinition |22| , |23 



9ab 

dip 



lab 



±^(^ + 3/2 



(2a) 
(2b) 
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the action becomes (in the so-called Einstein frame) 



S*= V=g 



# -5 <P*<P,b-2—a- + 2— 



c/ 4 x, (3) 



which is the action for general relativity (GR) containing a scalar field tp with 
the potential 



Our aim here is to exploit the results in previous work [17] to study 
the asymptotic properties of scalar-tensor theories of gravity with action ([!]) 
which under the transformations (0) transform to general relativity with a 
scalar field with the exponential potential given by 

V* = V e k(p , (5) 

where Vq and k are positive constants. That is, since we know the asymptotic 
behaviour of spatially homogeneous Bianchi models with action @ with the 
exponential potential ([|), we can deduce the asymptotic properties of the 
corresponding scalar-tensor theories under the transformations (||])[] (so long 
as the transformations are not singular!). In particular, we are concerned 
with the possible isotropization and inflation of such scalar-tensor theories. 

The outline of the paper is as follows. In section 2, we review the frame- 
work within which GR and a scalar field with a potential (Einstein frame) is 
formally equivalent to a scalar-tensor theory with a potential (Jordan frame), 
concentrating on both the exact and approximate forms for the parameters 
V and to in the Jordan frame. In particular, we discuss the explicit example 
of the Brans-Dicke theory with a power-law potential and we also discuss the 
conditions which lead to appropriate late-time behaviour as dictated by so- 
lar system and cosmological tests. In section 3, we then apply the conformal 
transformations to Bianchi models studied in the Einstein frame to produce 
exact solutions which represent the asymptotic behaviour of more general 
spatially-homogeneous models in the Jordan frame (for u = u , a constant). 
These Brans-Dicke models are self-similar and the corresponding homothetic 
vectors are also exhibited. We conclude with a discussion in section 4. 



lr The possible isotropization of spatially homogeneous scalar-tensor theories which get 
transformed to a model with an effective potential which passes through the origin and is 
concave up may be deduced from the results of Heusler EJ] 
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2 Analysis 




For scalar field Bianchi models the conformal factor in (|2a|) is a function of 
t only (i.e., = (f>(t)), and hence under (non-singular) transformations (0) 
the Bianchi type of the underlying models does not change (i.e., the metrics 
g a b and g* b admit three space-like Killing vectors acting transitively with 
the same group structure). In general, in the class of scalar-tensor theories 
represented by ([!]) there are two arbitrary (coupling) functions u(4>) and 
V((f>). The models which transform under (0) to an exponential potential 
model, in which the two arbitrary functions uj and V are constrained by 
and (f|), viz., 

J, l~o 

(6) 

is a special subclass with essentially one arbitrary function. Although only 
a subclass of models obey this constraint, this subclass is no less general 
than massless scalar field models (V = 0; see, for example fi~6f ) or Brans- 
Dicke models with a potential (oj = uq, constant), which are often studied 
in the literature. Indeed, the asymptotic analysis in this paper is valid not 
only for "exact" exponential models, but also for scalar-tensor models which 
transform under (|2]) to a model in which the effective potential is a linear 
combination of terms involving exponentials in which the dominant term 
asymptotically is a leading exponential term; hence the analysis here is rather 
more general (we shall return to this in the next section). For the remainder 
of this paper we shall not explicitly consider ordinary matter; i.e., we shall 
set the matter Lagrangians in (p]) and (|3|) to zero. Matter can be included 



in a straightforward way [24, 16, 25 . 



2.1 Exact Exponential Potential Models 

Scalar-tensor models which transform under (Q) to a model with an exact 
exponential potential satisfy equations ( pT5|) and (|J) with (|5]), viz., 



dtp yM0) + 3/2 

d(j) (j> U 



V«» = ^P- (8) 
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So long as the transformations (0) remain non-singular we can determine 
the asymptotic properties of the underlying scalar-tensor theories from the 
asymptotic properties of the exact exponential potential model. These prop- 
erties were studied in [|7j . We recall that the asymptotic behaviour depends 



crucially on the parameter k (in (|5D) which will be related to the various 
physical parameters in the scalar-tensor theory ([!]). 

In particular, in |l7j it was shown that all scalar field Bianchi models 
with an exponential potential ([5]) (except a subclass of the Bianchi type IX 
models which recollapse) isotropize to the future if k 2 < 2 and, furthermore, 
inflate if k 2 < 2; if k = these models inflate towards the de Sitter solution 
and in all other cases they experience power-law inflationary behaviour. If 
k 2 > 2, then the models cannot inflate, and can only isotropize to the future 
if the Bianchi model is of type /, V, VII, or IX. Those models that do not 
isotropize typically asymptote towards a Feinstein-Ibanez anisotropic model 
Tf\ . Bianchi Vllh models with k 2 > 2 can indeed isotropize |T7 but do 
not inflate, while generically the ever-expanding Bianchi IX models do not 
isotropize [p6|. 



Therefore, at late times and for each specific choice of uj(4>) both the 
asymptotic behaviour of the models and the character of the conformal trans- 
formation (0) may be determined by the behaviour of the scalar field ip at 
the equilibrium points of the system in the Einstein frame. Recently this 



behaviour has been thoroughly investigated |Ty|]. We shall summarize only 
those aspects relevant to our study. The existence of GR as an asymptotic 
limit at late times is also determined by the asymptotic behaviour of the 
scalar field; we shall return to this issue in section |2~3| . 

For spatially homogeneous space-times the scalar field <p is formally equiv- 
alent to a perfect fluid, and so expansion-normalized variables can be used to 
study the asymptotic behaviour of Bianchi models [IT], . The scalar field 
variable, is defined by 

* - 4", (9) 

where 6* is the expansion of the timelike congruences orthogonal to the sur- 
faces of homogeneity^. At the finite equilibrium points of the reduced system 
of autonomous ordinary differential equations, where \1/ is a finite constant, 
it has been shown |2S| that 9* = 9^/t*, where t* is the time defined in the 



2 Note that 0* > for all Bianchi models except those of type IX. 
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Einstein frame: 

dt*=±yfidt. (10) 

From equation (|9|) it follows that oc 1/t*, whence upon substitution into 
the Klein-Gordon equation 

dV* 

£ + + — = o, (ii) 

we find that at the finite equilibrium points 

<p(t*) =^ -|lnt*; k^O, (12) 

where (po is a constant. Hence, from equation fl2rj|) we can obtain as a 
function of £*, provided a particular cj(0) is given. From equation ( 110 ) we 
can then find the relationship between t* and t, and consequently obtain </> 
as a function of t, and hence determine the asymptotic behaviour of 4>(t) for 
a given theory with specific oo((f)) (in the Jordan frame). Specifically, we can 
determine the possible isotropization and inflation of a given scalar-tensor 
theory in a very straightforward way. 

As mentioned above, the behaviour determined from the key equation 
( |T2| ) is not necessarily valid for all Bianchi models. For the Bianchi models 
in which the phase space is compact, the equilibrium points represent models 
that do have the behaviour described by (|T^), as do the finite equilibrium 
points in Bianchi models with non-compact phase spaces. It is possible that 
the infinite equilibrium points in these non-compact phase spaces also share 
this behaviour, although this has not been proven. Finally, from equations 
(0) we note that since the asymptotic behaviour is governed by (|i~2[) , the cor- 
responding transformations are non-singular and this technique for studying 
the asymptotic properties of spatially homogeneous scalar-tensor theories is 
valid. 



2.2 An Example 

Suppose we consider a Brans-Dicke theory with a power-law potential, viz., 

u{4>) = u (13) 
V = (3(j) a (14) 
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(where (3 and a are positive constants), then (|2"T5|) integrates to yield 



?o exp 



where 

and hence (E) yields 



w = ±\/l}q + 3/2, 



where the critical parameter k is given by 



k 



a 



CO 



(15) 

(16) 
(17) 

(18) 



From [17| we can now determine the asymptotic behaviour of the models in 



the Einstein frame, as discussed in section |2.1| , for a given model with specific 
values of a and Q (and hence a particular value for k). 

The possible isotropization of the given scalar-tensor theory can now be 
obtained directly (essentially by reading off from the proceeding results - see 
subsection 3.1). For example, the inflationary behaviour of the theory can 
be determined from equations fl2a|), (|I0|) and (jTj). Let us further discuss the 
asymptotic behaviour of the corresponding scalar-tensor theories (in the Jor- 
dan frame). From equations (|T0D, (|TJ) and flTS] ) we have that asymptotically 



±(«-«.)(i + i5 



1 -2/(l+to) 



(19) 



where the ± sign is determined from (|10l). Both this sign and the signs of 
uj and 1 + kuj are crucial in determining the relationship between t* and t; 
i.e., as t* — > oo either t ±oo or t — > t and hence either — > or — > oo, 
respectively, as (y9 — >• — oo. 

2.2.1 A Generalization 



Suppose again that to = u;o, so that ( JT5J ) also follows, but now V is a sum of 
power-law terms of the form 



(20) 



n=0 
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where m > 1 is a positive integer. Then (|]) becomes 

m 

v* = J2p n r n ~ 2 

n=0 

= V p n exp(k n ip); k n = -^- — . (21) 

For example, if 

V = V + l -m^ 2 + A0 4 , 

then 

= y e" 2 ^ + -fh + Xe 2v/G) 

(with obvious definitions for the new constants), which is a linear sum of 
exponential potentials. Asymptotically one of these potentials will dominate 
(e.g., as if —> +oo, V* — > \e 2ip ^) and the asymptotic properties can be 
deduced as in the previous section. 



2.2.2 Approximate Forms 

In the last subsection we commented upon the asymptotic properties of a 



scalar-tensor theory with the forms for u and V given by ( [13|) and (0). Let 
us now consider a scalar-tensor theory with forms for u and V which are 
approximately given by ([X3| ) and ( fT^D (asymptotically in some well-defined 
sense) in order to discuss whether both theories will have the same asymp- 
totic properties. In doing so, we hope to determine whether the techniques 
discussed in this paper have a broader applicability. 

We assume that us and V are analytic at the asymptotic values of the 
scalar field in the Jordan frame in an attempt to determine whether their 
values correspond to the appropriate forms for ip and V* in the Einstein 
frame, namely whether (f — ► — oo and the leading term in V* is of the form 



Consider an analytic expansion for about = 0: 



oo 

u = (22) 

n=0 

oo 

V = (23) 

n=0 
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where all uj n and V n are constants. Using (^) we find, up to leading order in 
0, that for u ^ — 3/2 

<p — (po ~ cDln0, (24) 

so that </? — >■ ±00 (depending on the sign in (|T6"D) for — > 0. The potential 
in the Einstein frame is (to leading order) 

V^expj-^-^}. (25) 

Hence, the parameter k of ([!]) is defined here as = —2/uj. For cuq = —3/2 
we have 

(v?-v?o) 2 « 4c^0 (26) 

V - ? 16 "\ 4 , (27) 
(^-^o) 4 ' 1 ' 

so that ip -/-> — 00 as — > 0. 

Next, let us consider an expansion in 1/0, valid for — > 00: 

w = V^, (28) 
^ = ( 29 ) 



n=0 



For c^o 7^ — 3/2, the results are similar to the = expansion: 

(p — <f « — u;ln0 (30) 
V* - expj^^H, (31) 



where now 99 — ► =)=oo as — > 00. When cjq = — 3/2, we obt 



am 



(v^-^o) 2 - ^ (32) 

^ 16w? ' ( j 

It is apparent that the sign of u is important in determining whether — > 00 
or — > in order to obtain the appropriate form for (p, as was exemplified 
at the end of section 2.2. 
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Finally, in the event that ui and V are analytic about some finite value 
of 0, namely </>q, it can be shown that <p — > (fa as — > 0o- Hence, if one 
insists that to remain analytic as uj — > co>o in the limit of ^ — ► —00, then 
must either vanish or diverge, and the GR limit is not obtained. This would 
then suggest that if one imposed ip — > —00 for — > O then u would not be 
analytic about = 0q. 



2.3 Constraints on Possible Late-Time Behaviour 

In this paper we are concerned with the possible asymptotic behaviour of 
cosmological models in scalar-tensor theories of gravity. However, there are 
physical constraints on acceptable late-time behaviour (as t* — > 00; see equa- 
tion (|ToD ). For example, such theories ought to have GR as an asymptotic 
limit at late times (e.g., uj —>■ 00 and — > 0o) in order for the theories to 
concur with observations such as solar system tests. In addition, cosmo- 
logical models must 'isotropize' in order to be in accord with cosmological 
observations. 

Nordtvedt |7| has shown that for scalar-tensor theories with no potential, 
uj((p) — > 00 and u~ 3 duj/d(j) — > as t — > 00 in order for GR to be obtained in 
the weak-field limit. Similar requirements for general scalar-tensor theories 
with a non-zero potential are not known, and as will be demonstrated from 
the consideration of two particular examples found in the literature, not all 
theories will have a GR limit. 

The first example is the Brans-Dicke theory (u = ujq = constant) with a 
power-law self-interacting potential given by ( [14]) studied earlier in subsection 
2.2. In this case, is given by equations (|1^) and fllSD and the potential is 
given by ([14]), viz., 



V(<j>) = f3cj) a ; a = 2 t k Jlu + 3/2. 



The a = 1 case for FRW metrics was studied by Kolitch p9| and the a = 2 
(k = 0) case, corresponding to a cosmological constant in the Einstein-frame, 
was considered for FRW metrics by Santos and Gregory |30|. Earlier we 
considered whether anisotropic models in Brans-Dicke theory with a potential 
given by (|I4]) will isotropize. Assuming a large value for u> , as suggested by 



solar system experiments, we conclude that for a wide range of values for 
a the models isotropize. However, in the low-energy limit of string theory 
where uq = — 1 the models are only guaranteed to isotropize for 1 < a < 3. 
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Substituting (|TJ) in fll^) we get 



^' r » ±M ^WaTW (34) 

Now, substituting the above expression into equation ([5]), we obtain £* as a 
function of t and hence we obtain 

(35) 

Depending on the sign, we deduce from this expression that for large t the 
scalar field tends either to zero or to infinity and so this theory, with the 
potential given by ([14]), does not have a GR limit. 
In the second example we assume that 

"(0 + | = 7T^2> (36) 
2 {(j) - O ) 

where A is an arbitrary positive constant. This form for ou((f)) was first 
considered by Mimoso and Wands Jl6| (in a theory without a potential). 
Now, we obtain 

= o + fi e T7i, (37) 

where B is a constant, and the potential, defined by equation (H), is given 
by 

y(0) = yo0 2 (0-0 o r v/3fc . (38) 

As before, at the equilibrium points we can express as a function of £*, 
which then allows us to compute i as a function of t*. At late times we find 
that 

0~0o+^, (39) 

where (3 is a constant whose sign depends on k, u and the choice of one 
of the signs in the theory. What is important here is that in this case, at 
late times, we find that the scalar field tends to a constant value for (3 < 0, 
thereby yielding a GR limit. In both of the examples considered above, the 
conformal transformation for the equilibrium points is regular. 

Of course, these are not the only possible forms for a variable o;(0). For 
example, Barrow and Mimoso [^] studied models with 2u((j>) + 3 oc <f) a (a > 
0) satisfying the GR limit asymptotically. (The GR limit is only obtained 
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asymptotically as — > oo, although for a finite but large value of the theory 
can have a limit which is as close to GR as is required). However, by studying 
the evolution of the gravitational "constant" G from the full Einstein field 
equations (i.e., not just the weak-field approximation), Nordtvedt |7|, |31| has 
shown that 

G /3 + 2oA A 2sJ 
G ~ ~ 



where uj' - 
and cu'u~ 3 



duj/d(f) 
- 0). 



4 + 2a; 7 V" ' (3 + 2a;) 2 J 

so that the correct GR limit is only obtained as uj — ► oo 
Torres [32] showed that when 2o>(0) + 3 oc a , G(t) 
decreases logarithmically and hence G — > asymptotically. In the above 
work, no potential was included. For a theory with 2a; (0) + 3 oc Q and with 
a non-zero potential satisfying equation (|6]) we have that 

(a 7^ 0; A and B constants), so that 

V{(f) = V (j) A e B ^ /a . 



A potential of this form was considered by Barrow [33| . 

Finally, Barrow and Parsons [^3] have studied three parameterized classes 
of models for oj{4>) which permit uj — ► oo as — > <po (where the constant 
0o can be taken as evaluated at the present time) and hence have an 
appropriate GR limit; 



2o;(0) + 3 = 2B\ |1 



{%%) 2cj(0) + 3 = 5 2 2 |/n(0/0 o )| 



-2\f\ 



[III 



2cu(0) + 3 = B\ 



1 - 



( 1 ^ 

(«>5), 

(V/3). 



Other possible forms for u(4>) were discussed in Barrow and Carr and, 
in particular, they considered models (i) above but allowed a < in order 
for a possible GR limit to be obtained also as — ► 00. Schwinger [36| has 
suggested the form 2a; (0) + 3 = B 2 /(f) based on physical considerations. 
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3 Applications 



Let us exploit the formal equivalence of the class of scalar-tensor theories ([I]) 
with u)(<f>) and V((f>) given by 

u(<P) = u Q , V{<f>) = P<f> a , (40) 

with that of GR containing a scalar field and an exponential potential @. 
Indeed, since the conformal transformation fl2a]) is well-defined in all cases 
of interest, the Bianchi type is invariant under the transformation and we 
can deduce the asymptotic properties of the scalar-tensor theories from the 
corresponding behaviour in the Einstein frame. Also, we have that 

k = ^— , iu 2 = cu + -. 41 
UJ 2 

We recall that at the finite equilibrium points in the Einstein frame we 
have that 

6* = ex 1 , (42) 

<p{t*) = ^ -^ln(r), (43) 

where 

9* = 1 + —e kvo . (44) 

Integrating equation ( ppj ) we obtain 

<f>(f) = rfexp (w-V(**)) = (PoK 2 ^, (45) 

where the constant 0o = dexp(ip Q /u) and we recall that t and t* are related 
by equation (0), and equation ( pa]) can be written as 

^ = (46) 

3.1 Examples 

1) All initially expanding scalar field Bianchi models with an exponential 
potential @ with < k 2 < 2 within general relativity (except for a sub- 
class of models of type IX) isotropize to the future towards the power-law 
inflationary flat FRW model |25], whose metric is given by 



ds 2 = -dtl + t! k2 (dx 2 + dy 2 + dz 2 ) . (47) 
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In the scalar-tensor theory (in the Jordan frame), is given by equation (|5) 
and from (Hot) we have that 



ds 2 ST = {ds 2 } . (48) 

Defining a new time coordinate by 

l + kuj k'UJ — 1 

T = cU^; c=__0 o 5 (49) 
1 + kuj 

(where kuo + 1 ^ 0; i.e., a ^ 1), we obtain after a constant rescaling of the 
spatial coordinates 

ds 2 ST = -dT 2 + T 2K (dX 2 + dY 2 + dZ 2 ) , (50) 

where 

k 2 + 2koo 
= k 2 (l + ku)' 
Finally, the scalar field is given by 

2 — 2 — 2 

= oC T+fe^Ti+^ = O T~. (51) 

Therefore, all initally-expanding spatially-homogeneous models in scalar- 
tensor theories obeying (E0|) with < (a — 2) 2 < 2u + 3 (except for a 
subclass of Bianchi IX models which recollapse) will asymptote towards the 
exact power- law flat FRW model given by equations (|5(J) and fl5"ID, which 



will always be inflationary since K = 7~]rfcz%) > ^ [ n °t e that whenever 
2cj > (a - 2) 2 - 3 = a 2 - 4a + 1, we have that 1 + a + 2u > a 2 - 3a + 2 = 
(a-l)(a-2)]. 

When k 2 > 2, the models in the Einstein frame cannot inflate and may 
or may not isotropize. Let us consider two examples. 

2) Scalar field models of Bianchi type VI^ with an exponential potential 
(H) with k 2 > 2 asymptote to the future towards the anisotropic Feinstein- 
Ibanez model |27J given by (m ^ 1) 

ds 2 = -dtl + a 2 (tl P1 dx 2 + t 2p2 e 2mx dy 2 + t 2 / 3 e 2x dz 2 ) , (52) 
where the constants obey 



Pi = 


1- 










P2 = 


k 2 








P3 = 


k 2 



(A; 2 - 2)(m 2 + m)' 
2(m 2 + 1) / 
(£; 2 -2)(m + l)\ 
2(m 2 + 1) ) 



(53) 
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In the scalar-tensor theory (in the Jordan frame), is given by equation 
(|45p and the metric is given by (|48|). After defining the new time coordinate 
given by (|49|), we obtain 

ds 2 ST = -dT 2 + Al (T 2qi dX 2 + T 2q2 e 2mX dY 2 + T 2qs e 2X dZ 2 ) , (54) 

where 

( < = 1 ' 2 ' 3 ); Al = al4>- Q l c^\ (55) 

and Y and Z are obtained by a simple constant rescaling (and X = x). 
Finally, the scalar field is given by equation ([51]). 

The corresponding exact Bianchi VI^ scalar-tensor theory solution is 
therefore given by equations ([51]) and (|54"D in the coordinates (T, X, Y, Z). 
Consequently, all Bianchi type VI^ models in the scalar-tensor theory satis- 
fying equations fl4"0[) with (a — 2) 2 > 2u> + 3 asymptote towards the exact 
anisotropic solution given by equations flsT]) and (|5"4[) . 

3) An open set of scalar field models of Bianchi type VII^ with an 
exponential potential with k 2 > 2 asymptote towards the isotropic (but non- 
inflationary) negative-curvature FRW model [|I7[] with metric 

ds 2 = -dtl + tlda 2 , (56) 

where da 2 is the three-metric of a space of constant negative curvature. 
Again, is given by (^) and the metric is given by (0), which becomes 
after the time recoordinatization 



ds 2 ST = -dT 2 + C 2 T 2 da 2 , (57) 

r - n 2 

where C 2 = 4>q 1 c~ 2 = . This negatively-curved FRW metric is equiv- 

alent to that given by (0). Finally, the scalar field is given by equation 
&■ 

Therefore, when (a — 2) 2 > 2^0 + 3, there is an open set of (BVII^) scalar- 
tensor theory solutions satisfying equations ( ^0|) which asymptote towards 
the exact isotropic solution given by equations (|5l|) and fl57l) . 

Equations (|4"3"P and (H3) and the resulting analysis are only valid for 



scalar-tensor theories satisfying (f40|). However, the asymptotic analysis will 
also apply to generalized theories of the forms discussed in subsections 2.2.1 
and 2.2.2. Finally, a similar analysis can be applied in Brans-Dicke theory 
with V = |B7H. 
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3.2 Self-Similarity 

All three attracting scalar-tensor theory solutions in the last subsection are 
self-similar; metric (50) admits the homothetic vector (HV) X = T-Jj^ + 



(1 - K) {X-^ + Y-^p + Z-gz}, metric (|§ admits the HV X = T ^ + (1 - 

q 2 )Y-7^: + (1 — q 3 )Z-^, and metric fl57j) admits the HV X = T-^. Of course, 
all three solutions in the corresponding general relativistic model (i.e., in the 
Einstein frame) are self-similar. Let us show that this is always the case; i.e., 
all scalar-tensor solutions obtained in this way are self-similar. 

In |L7j it was shown that the cosmological solutions corresponding to the 
finite equilibrium points of the "reduced dynamical system" of the spatially 
homogeneous scalar field models with an exponential potential are all self- 
similar. Let g* ah be the metric of such a solution and X* the corresponding 
HV; hence we have that 

LXJL = 2 ^ (58) 
where L denotes the Lie derivative along X*. In the coordinates in which 
0* = from £x„#* = -9* we find that [gTj 



X, = t .A + ^(* Z )A (59) 

Now, the metric g ab in the corresponding scalar-tensor theory is given by 
), where the scalar field is given by ( ^5|) , viz., 

<KO = o t; 2/to (60) 

(or by (|5lD in terms of the time coordinate T). We emphasize that this power- 
law form for is only valid for scalar-tensor theories that obey conditions 
PD|). Hence, from equations (|58"])-(|6"0|) it follows that 



c XJab = x^r^g^+r'CxXb 

= f J? (V^) 9* ab + 2<f>o 1 tPg* ab 



kuo 

= 2cg ab , (61) 

where the constant c is given by c = (1 + kuj)/{kuj). That is, X = X* is 
a homothetic vector for the spacetime with metric g ab and consequently the 
corresponding scalar-tensor theory solution is self-similar. 
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3.3 The Special Case a = 1 

In the analysis above we have omitted the special case a = 1 (i.e., ku = 
— 1). This case is degenerate as we will now demonstrate. Let the general 
relativistic metric be defined by 

ds 2 = -dt\ + y^dx^dx" . (62) 

First, suppose we take ku = — 1 in ([PJ) and define a new time coordinate by 

T = 0~Mn(t«), (63) 

then the metric (^) becomes 

ds 2 ST = -dT 2 + 0o 1 exp (^JfoT^y^dx^dx", (64) 



where <f){T) = 0oexp (2^/<j)oT^. Now, from equation (^) we obtain 

C>X*9ab = 0; (65) 



i.e., in this case X = X* is a Killing vector (KV) of the spacetime (|64]). Since 
the KV X is timelike, the spatially homogeneous metric (|64D admits four KV 
acting simply transitively and hence the resulting spacetime is (totally - i.e., 
four-dimensionally) homogeneous. 



All known non-flat homogeneous spacetimes are given in table 10.1 in |38 
hence the metric (^4|) is given by one of those spacetimes in this table repre- 
senting an orthogonal spatially homogeneous metric with a diagonal Einstein 
tensor (representing a perfect fluid spacetime or an Einstein spacetime with 
a cosmological constant) - all of these metrics are indeed known f38| . In the 



case when metric ( |B"4] ) is the flat Minkowski metric, the corresponding general 
relativistic spacetime (|62|) is de Sitter spacetime. However, this corresponds 
to the degenerate case in which 

6* = 9q, a constant; 

this is the only possibility in which equation (^) is not valid and hence the 
resulting analysis does not follow. This degenerate case corresponds to k = 
in (H) (i.e., V* = Vo, a constant); since kuj = —1 this corresponds to uj — > oo 
or u)q —>■ oo (in which case GR is recovered from the scalar-tensor theory 
under consideration). 
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Finally, if a = 1 in ( fH| ) (i.e., V = then the action ([I]) becomes 



S 



<P{R-2(3) 



+ 2£ r 



d x, 



which is equivalent to that for Brans-Dicke theory incorporating an additional 
constant (3. Under the conformal transformation and field redefinition (Q) the 
action becomes that for general relativity with a cosmological constant (and 
additional matter fields), and from the cosmic no-hair theorem it follows that 
all spatially homogeneous models (except for a subclass of Bianchi type IX) 
asymptote to the future towards the de Sitter model p9|. 



4 Conclusions 

In this paper we have studied the asymptotic behaviour of a special subclass 
of spatially homogeneous cosmological models in scalar-tensor theories which 
are conformally equivalent to general relativistic Bianchi models containing 
a scalar field with an exponential potential by exploiting results found in 
previous work |T7|| . 

We illustrated the method by studying the particular example of Brans- 
Dicke theory with a power-law potential and various generalizations thereof, 
paying particular attention to the possible isotropization and inflation of 
such models. In addition, we discussed physical constraints on possible late- 
time behaviour and, in particular, whether the scalar-tensor theories under 
consideration have a general relativistic limit at late times. 

In particular, several exact scalar-tensor theory cosmological models (both 
inflationary and non- inflationary, isotropic and anisotropic) which act as at- 
tractors were discussed, and all such exact scalar-tensor solutions were shown 
to be self-similar. 

This is related to the previous work of several authors. Specifically, Chau- 
vet and Cervantes-Cota studied isotropization in Brans-Dicke gravity 



including a perfect fluid with p = (7 — 1)^. They examined whether the 
anisotropic models contain an FRW model as an asymptotic limit, which is 
how they defined isotropization. For Bianchi models of types J, V and IX, 
they found exact solutions in these cosmologies which can isotropize to the 
future, depending on the values of 7 and uj and two other arbitrary constants. 
Furthermore, Mimoso and Wands |16| also studied scalar-tensor models with 



variable u without a self-interacting potential V but coupled to barotropic 
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matter. Regarding the possible isotropization of the cosmological models 
(meaning here that the shear of the fluid becomes negligible), they concen- 
trated on models of Bianchi type / and first discussed constraints on a fixed 
uj = co>o model necessary for isotropization at late times. In the particular 
case of a false vacuum (p = — //), they showed that the de Sitter solution 
is the late-time attractor of the model. They then proceeded to examine 
arbitrary uj((p) Bianchi type / cosmologies and showed that if a solution is 
to asymptote towards a GR limit (i.e., uo — > oo) then it must also isotropize. 
Their paper also discussed initial singularities in models of other Bianchi 
types. 

The work in this paper can be generalized in a number of ways. In 
particular, more general scalar-tensor theories can be considered and more 
general (than spatially homogeneous) geometries can be studied. For exam- 
ple, the more general class of inhomogeneous G2 models could be considered 



Mj |4], [|!| in which there exists two commuting spacelike Killing vectors. 
The motivation for studying G 2 cosmologies is that there is some evidence 
that the class of self-similar G 2 models plays an important role in describ- 
ing the asymptotic behaviour of more generic general relativistic scalar field 
models with an exponential potential [cf . 42] ; in this way, we may be able to 
find special scalar-tensor G2 cosmological models that describe the asymp- 
totic properties of more general scalar-tensor cosmologies. Some potential 
problems that exist in this more general context is that since 0, and hence 
the transformation (|2"a|), depends on both time and one space variable, the 
transformation will be singular (at least for certain values of the space 
variable) and the classification of G2 models may not be preserved under 
such a transformation. 
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